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Abstract 

We study the quantum entanglement and quantum phase transition (QPT) of the anisotropic 

O ■ spin-1/2 XY model with staggered Dzyaloshinskii-Moriya (DM) interaction by means of quantum 

QJ . 

^ , renormalization group method. The scaling of coupling constants and the critical points of the 

' system are obtained. It is found that when the number of renormalization group iterations tends 

> ' 
c/2 ' 

to infinity, the system exhibit a QPT between the spin-fluid and Neel phases which corresponds 
with two saturated values of the concurrence for a given value of the strength of DM interaction. 



' The DM interaction can enhance the entanglement and influence the QPT of the system. To gain 

O ' further insight, the first derivative of the entanglement exhibit a nonanalytic behavior at the critical 



point and it directly associates with the divergence of the correlation length. This shows that the 
correlation length exponent is closely related to the critical exponent, i.e., the scaling behaviors of 



O ' the system. 

O 
O 
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I. INTRODUCTION 



In the quantum systems, entanglement is a pure quantum correlation, which is the fun- 
damental difference between quantum and classical physics In recent years, quantum 
entanglement has attracted much attention in quantum information theory because of its 
importance in developing the idea of quantum computers and other quantum information 
devices 3|. It has also been realized as a crucial resource to process and send information 
in different ways, such as quantum teleportation, quantum cryptography, and algorithms 
for quantum computations |#^6|. In the condensed-matter physics, it is very significant to 
discuss the relation between entanglement and quantum phase transition (QPT) which has 
been attracted many researchers to inverstigate 7H9|. 

For investigating the properties of many-body systems, the renormalization-group (RG) 
method is applied. In the past several decades, much effort had been investigated in many 
spin systems using this method. Real space renormalization group method was applied 
;o discuss the critical points and phase diagrams of Heisenberg and Blume-Capel models 
10-13|. These properties of some models were also discussed by Monte Carlo RG [l^ Il5 |. 



The density-matrix RG method which is a powerful numerical method is used to study 
the ground and low-lying states properties of low- dimensional lattice models. It has been 
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fed successfully to lots of strongly correlated systems in ID as well as 2D systems 
19|. Recently, the pairwise entanglement of the system is studied by the quantum 



renormalization-group (QGR) method which plays an important role in QPT 20|, |21|. The 



spin-1/2 Ising and Heisenberg models were investigated by the same method and it is found 
that the systems exist QPT and nonanalytic behavior, such as the discontinuity in the 



quantum critical points 



22 



25| . For getting the accurate results, the XXZ model with next- 



nearest-neighbor interactions are investigated |26l . |27[ . It is found that the tri-critical point 
and the phase diagram of the sysytem can be obtained. 

Some spin models can be supplemented with a magnetic term which is called 
Dzyaloshinskii-Moriya (DM) interaction arising from the spin-orbit coupling. The DM inter- 



action, which was first proposed by Dzyaloshinskii and Moriya about half century ago|28ll29|. 
can influence the phase transition and the critical properties of some systems. The relevance 
of antisymmetric superexchange interaction which describes quantum antiferromagnetic sys- 
tem was introduced by Dzyaloshinskii. Moriya found that such interaction arises naturally 



2 



in the perturbation theory in magnetic systems with low symmetry. The form of DM interac- 
tion for two spins 5^ and Sj is X S^^ . Ising and XXZ models with DM interaction 



were studied in Ref. 



22 



23|. The results are that the critical points of the systems are 



obtained and divided the systems into two phases, i.e., spin-fluid and Neel phases. At the 
critical point, the nonanalytic behavior of the first derivative of the entanglement and the 
scaling behavior of the systems are also gotten. 

The quantum entanglement and QPT of the spin-1/2 XY model with staggered DM 
interaction are discussed by using QRG method. We find that the stable and unstable fixed 
points of the system can be obtained and the phase transition point changes as the DM 
interaction increases.. The concurrence is calculated which is influenced by the anisotropy 
parameter and DM interaction. The concurrence trends two fixed values which associate 
with the phases of the system as the number of RG iteration increases. Furthermore, the 
first derivative of the concurrence shows nonanalytic behavior at the critical point which 
has relation with the correlation length. This paper is organized as follows. In Sec. [Tll we 
apply QRG method to investigate the model and obtain the fixed points. The concurrence 
is introduced to measure the entanglement and we analysis it in order to get more insights 
about the critical features of the model in Sec. Illli We summarize in Sec. IIV[ 



II. QUANTUM RENORMALIZATION GROUP OF THE MODEL 

The mode elimination or the thinning of the degrees of freedom followed by an iteration, 
which reduces the number of lattices step by step until reaching a more tractable circum- 
stance, is the main idea of RG method. RG is a proper method to give the universal behavior 
at long wavelengths, it includes many methods, such as decimation, bond-moving and cu- 
mulant expansion. In this paper, the Kadanoff's block approach is implemented where we 
have consider three sites as a block (marking as 1-2-3). Generally speaking, this method 
includes three steps. Firstly, the system is divided into blocks and the Hamiltonian of each 
block can be exactly diagonalized and solved. Then, the projection operator is builded by 
the lower eigenvectors. And finally, the full Hamiltonian is projected onto these eigenvectors 
to obtain the effective Hamiltonian which acts on the renormalized subspace, i.e., the RG 



equations 



30 



3l|. 



The Hamiltonian of ID anisotropic XY model with staggered DM interaction on a peri- 
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odic chain of sites can be written as 

N 



J 

where J is the nearest exchange couphng constant, 7 is the anisotropy parameter, D is the 
strength of DM interaction in the direction of z, and {a = x, y) are Pauh operators of 
the ith site. For the Hq, when 7 = D = 0, the model becomes isotropic XX model; when 
7 = 1 and D = 0, it turns into the Ising model which was exactly solved in Ref. j32| . 

The initial Hamiltonian Hq acts on the effective Hilbert space and then the effective 
Hamiltonian H'^^ can be gotten. The essential criterion for RG is that H"^ have similar 



structure to Hq, but we can not get this, i.e., the signs of o"fcrf+i and cTiCrf^i terms are 



changed. To avoid this and produce a self-similar Hamiltonian, we imp 
around the x axis for all even sites and leave all odd sites unchanged 
transformed Hamiltonian is obtained as follows, 

N 

4' 



+1 

ement a vr rotation 



25| . Therefore, the 



i=l 



For the Kadanoff 's block approach, H can be written as 

H = H'' + i/^^, (3) 

where H^ is the block Hamiltonian and H^^ is the interblock Hamiltonian. The explicit 
forms of H^ and H^^ are 

N/3 

H' = Y.hf. (4) 

1=1 

N/3 

= I E [(1 + ^) <3^r+i,i - (1 - 7) <3<i,i + D «3<,,, + <3^r+i,i)] ^ (5) 
1=1 

where the Ith block Hamiltonian is 

hf = ^ [(1 + 7) «1<2 + <2<3) - (1 - 7) + <2<3) 

+D {afX, + alal, + af^.a^ + a^al,)] . (6) 



In terms of matrix product states [33!], we can get the eigenvalues and eigenvectors by 
exactly solving hf . The ground states which are doubly-degeneracy are useful to construct 
the projection operator and calculate the entanglement in later. In the standard basis 



{Ittt) , Itti) 7 Itit) > Itii) J litt) , \ltl) , lilt) > liii)}, and the degenerate ground states are 
given by 



l<^o) 



l<^o)' 



Itti) + 



7 



Di + 1 



Itit) + 



-V2{1 - Di) 



Ittt) + Itii) 



V2{Di + l) 
727 



litt) + liii) 



liti) + liit) 



, (7) 



(8) 



q q 

where q = ^J\ + D"^ + 7^, |'|") and |i) are the basis vectors of in itself representation. The 
energy corresponding to the ground states is 



eo 



-g/V2. 



We keep the ground states of hf to define the effective site. The effective Hamiltonian 
/f^ff and the Hamihonian H have in common the low lying spectrum [3^. An exactly 
implementation of this is given by the following equation, 



(9) 



Ar/3 



In the preceding equation, T = H ^0 is the projection operator of the system and the 
specific form of Tq is 



1=1 



^o = l^)^(v^o| + |^),(^o|' 



(10) 



where |f|^)j and |-l|); are the renamed states of each block in the effective space, which can be 
seen as a different spin-1/2 particle. Here, we consider only the first-order correction in the 
perturbation theory. By using the Eq. ([3]), the Eq. ([9]) can also be written as 



The Pauli matrices in x and y directions of renormalization are obtained as follows. 



(11) 



T},^ain=f^,a'r + u,ay (j = l,2,3). 



(12) 
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where 



?! — S3 — -^===^=^1 U — ^3 



7D 7 - 1 - ^2 

1^1= 1^3= , ,^ ^ '^l = i^3 



We substitute them into Eq. f ITT]) and obtain the effective Hamiltonian, 

, 7V/3 
A:=l 

where 

1 + ^2^37^ , 37 + 3^2^7^ , 
^ = 2? ^ = l + D2 + 3y > ^ = (15) 

Because the antisymmetric is the special property of DM interaction, i.e., Z^jj = — Z^j j, 
the stable and unstable fixed points can be gotten by solving 7' = 7. The stable fixed points 



locate at 7 = 00 and 7 = 0, the unstable fixed point is 7 = ±-\/l + -D^ which separates the 



spin-fluid phase, 7 = and 7 = 00, from the Neel phase, < I7I < vT+Z^. 



III. ENTANGLEMENT ANALYSIS 



There are many measures for pairwise entanglement 35l-l38l|. In this section, we would 
like to calculate the concurrence of pure state where we consider one of the degeneracy 
ground states. The density matrix of a ground state is composed, i.e., 

P = \fo) {fol ■ (16) 

The results of using [(/^q) to construct the density matrix will be same as Eq. flTB]) . 

There are two cases to define the concurrence for a three-site block, (z) The concurrence 
between sites 1 and 3 is obtained by summing over the degrees of freedom of the middle site 
2. (ii) We trace over the site 1 or 3 and get the concurrence between the middle site 2 and 
the other one. Without loss of generality, we consider the case (i). In the standard basis 
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{|tt) ) Iti) ) lit) 7 III)} J the reduced density matrix pi3 for sites 1 and 3 can be gotten in Eq. 

dsi), 



Pi3 = Tr2 [p] 



1 + 
2g2 



7^ 
1+D2 








h \ 

D+i 














2{1+D2) 





9^ 







2(1+D2) 


2(1+D2) 


— 47 
in 








1 / 



(17) 



Ci3 denotes the concurrence between the sites 1 and 3 which is defined as 



C 



13 



max {Ai — A2 — A3 — A4, 0} 



(18) 



where A^ {k = 1,2,3,4) are the sc 
order. pi3 is the spin-fiipped state 



uare roots of eigenvalues of = puPis in descending 
38| and its definition is 



Pi3 = {(^1 ® crl) pI^ {af (g) (t|) 



(19) 



where p*3 is the complex conjugate of P13. The value of C13 ranges from to 1, if C13 = 
or 1, the system is in an unentangled or a maximally entangled state, else it corresponds to 
a partial entangled state. The square eigenvalues of R are gotten. 



Ai — -, A2 — — , A3 — A4 — 0. 

2 q 



(20) 



Thus, the concurrence is obtained as follows, 



Ci3 = max{Ai - A2, 0} 



1 -fVl + D^ 



(21) 



2 q 

It is easy to see that C13 is infiuenced by 7 and D. For three-site model, we plot C13 
versus 7 for different values of D in Fig. 1. From the figure, it is found that the entanglement 
is a fixed value regardless of any value of D at 7 = or infinity. That is to say there are not 
phase transition for the XX model. Moreover, the DM interaction plays an important role 
and enhances the entanglement of the system when 7 is small comparing to D; the effect of 
anisotropy parameter for the entanglement is more important than DM interaction when 7 
is large. 

The purpose of QRG is that the full properties of the model enter a few sites through the 
renormalizing of coupling constants. The renormalization of the strength of DM interaction 
and anisotropic parameter are obtained which are contribution to the concurrence of the 
system. For a fixed value of = 1, the graph for C13 and 7 is plotted in Fig. 2. It 



reveals that as the number of QRG iterations increases, the concurrence develops three 
rather different features which are separated by 7 = and 7 = -\/2. When the number of 
RG iterations is very large, i.e., the system is infinity, the value of C13 is zero corresponding 
to Neel phase for 7 ranging from — -\/2 to except the point of 0; at 7 = 0, the system is 
the spin-fluid phase corresponding to the maximum value of the concurrence; for I7I > v^, 
the concurrence of the system slowly increases with I7I increasing and finally reaches to the 
maximum value. The system is in the same phase for 7 = and 7 = 00. For this, it can be 
explained from the Hamiltonian, i.e., the anisotropic XY model changes the isotropic XX 
model for 7 = and 7 00. 

Further insight, when the system is large enough, the first derivation of the entanglement 
shows the nonanalytic behavior at the critical point. Because dCis/d'-f is an even function 
of 7, we have only plotted dCi^/dj in 7 > for L> = 1 in Fig. 3. When the QRG iteration 
trends infinity, there are a minimum and a maximum values for each plot. From the diagram, 
it is found that the singular behavior of the concurrence becomes more pronounced at the 
thermodynamic limit. For the inset of diagram, there is a maximum value for each curved 
shape which verges to zero as the system becomes large. These also manifest that there are 
same properties for 7 = and 7 = 00. For a more detailed analysis, the positions of the 
minimum or the maximum of dCis/dj with the size of system increasing are given in Fig. 4. 
It can be seen that it shows a linear behavior. That is to say the position of the minimum or 
maximum of dCis/d^ tends toward the critical point 7 = vT+~D^. The critical exponent 9 
for this behavior is dC 13/ d^yl^^.^ or dCis/d^yl^^^ ~ xhese results justify that 9 is the 

reciprocal of the correlation length exponent v closing to the critical point, i.e., 9 = 1/v. 

Next, we investigate how the concurrence versus D changes for a fixed value of 7. The 
concurrence with D changing for different iterations is depicted in Fig. 5. It is clearly seen 
that the concurrence becomes large with D increasing and then until to 0.5 at last. But 
the larger the size of system is, the slower the increasing tendency of the concurrence is. 
This implies that the system can not occur QPT by changing the DM interaction. The first 
partial derivative of concurrence for D is clearly discussed in Fig. 6. There is not suddenly 
mutative in the diagram, but the maximum value of each plot becomes small with the size 
of the system largening. While we can also get the linear behavior between the maximum 
of dCi^/dD and the corresponding size of system which is plotted in Fig. 7. The exponent 
for the curve is dCis/dD\D^^^ ~ ^-o.98 fg^^g ^^^^ ^j^g Ising-hke universality class. The 
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position of the maximum of concurrence tends to infinity as the iteration of QRG increases, 
in other words, the stable fixed point — > oo is reached. 

IV. CONCLUSIONS 

The relation between the concurrence as a measure of quantum correlations and QPT was 
obtained. The anisotropy and DM interaction parameters determined the phase diagrams 
of the model. If the anisotropy parameter was small, the concurrence depended on the DM 
interaction, else the anisotropy parameter played an important role. As the number of RG 
iterations increased, the concurrence developed two different values which separated two 
phases for a given D, i.e., spin- fluid phase and Neel phase. For Neel phase, the larger the 
value of D was, the wider the width of value of 7 was. The critical behavior was described 
by the first derivative of the concurrence of the blocks. The scaling behavior characterizes 
how the critical point of the model was touched as the size of system increased. The critical 
exponent had relation with the correlation length exponent in the vicinity of the critical 
point. This implied that the quantum critical properties of the model were closely associated 
with the behavior of the entanglement. The concurrence increased slowly as the the size of 
the system became large, but the tendency of concurrence was unanimous, namely trending 
a fixed value. 
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Figure Captions 

Fig 1. (Color online) Concurrence between the first and third sites of a three-site model 
in terms of anisotropy for different values of the strength of DM interaction. 

Fig 2. (Color onhne) Representation of the evolution of the concurrence in terms of QRG 
iterations at a fixed value oi D — 1. The inset of diagram of 7 ranges from — \/2 to V2. 

Fig 3. (Color online) First derivative of concurrence and its manifestation toward diverg- 
ing as the number of QRG iterations increases (Fig. 2). The value of the inset of diagram 
manifestation toward zero. 

Fig 4. (Color online) The logarithm of the absolute value of minimum or maximum, 
ln{\dCis/ d'y\^rn\) , versus the logarithm of chain size, ln(A'"), which shows a scaling behavior. 
Each point corresponds to the minimum or maximum value of a single plot of Fig. 3. 

Fig 5. (Color onhne) Representation of the evolution of the concurrence in terms of QRG 
iterations at a fixed value of 7 = \/2- 

Fig 6. (Color online) First derivative of concurrence and its manifestation toward zero 
as the number of QRG iterations increases (Fig. 5). 

Fig 7. (Color onhne) The logarithm of the absolute value of maximum, 
ln{\dCis/ dD\D^gjx.\) : versus the logarithm of chain size, In(A^), which is hnear and shows 
a scaling behavior. Each point corresponds to the maximum value of a single plot of Fig. 5. 
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